Fig.l, dispersion(c=l,m=2,d=0.02,u=2) 
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We investigate the dynamical chiral symmetry breaking ( DCSB ) and superconductivity in a 
supersymmetric model at finite temperature and density. We employ the TV = 1 four-dimensional 
generalized supersymmetric Nambu— Jona-Lasinio model ( M = 1 generalized SNJL4 ) with a chem- 
ical potential as the model Lagrangian, and select the gauge freedom as U(l). In order to realize 
the DCSB and BCS-type superconductivity in this model, we introduce a SUSY soft mass term. 
Under the finite-temperature Matsubara formalism, the effective potential and the gap equations are 
derived in the flamework of the large- iV expansion. The finite-density effect in the DCSB is shown 

QQ ' by the critical coupling. The roles of both the boson and fermion sectors in the superconductivity 

f*"^ ' are examined by the quasiparticle excitation spectra and the gap equations. 

o 
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5h ! 

<^ ■ I. INTRODUCTION 

' Until now, the BCS ( Bardeen-Cooper-Schrieffer ) theory of superconductivity has influenced a large part of theoret- 
ical physics [1,2,3,4,5,6,7,8]. In the investigations of color superconductivity ( CSC ), the BCS method in relativistic 
field theory is widely used [6]. Because of its fundamental importance, it is interesting to construct a supersymmetric 
k> , theory of superconductivity. 

Recently, several attempts to examine the theories of CSC from the viewpoint of supersymmetry ( SUSY ) appeared 
in literature [9,10]. According as the results of the nonperturbativc method of SUSY gauge theories [11], Ref. [9] 
discussed the symmetry breaking patterns, especially the breaking of baryon density symmetry U(1)b of the J\f = 1 
SQCD ( supersymmetric QCD ) with a nonzero chemical potential fj, and a SUSY breaking mass A for squarks. 
, However, any gauge-symmetry-breaking two-body pairs like diquarks were not treated in Ref. [9]. Furthermore, 
because the validity of the exact results of SUSY gauge theories to the problem of CSC is not clear, the method 
of Ref. [9] could apply to the situation of fi < A <C A^sqcd, where [i and A ( both break SUSY explicitly ) could 
be regarded as a small perturbation to a SUSY gauge theory. Under a similar context, the authors of Ref. [10] 
proposed a toy model for giving diquark-like condensates. They introduced an SO(N) gauge interaction between 
Oh! quark superfields stronger than that of SU(3 C ), and then they argued the SO(N) gauge dynamics gives diquark-type 
condensations in their model. This SO(N) gauge interaction seems artificial. 

The purpose of this work is different from that of Refs. [9,10]. We investigate the dynamical chiral symmetry 
breaking ( DCSB ) and superconductivity in a SUSY condensed matter at fi ~ O(A). fi is the characteristic energy 
scale of a condensed matter system, while A gives the SUSY breaking scale. Roles of superpartners in the DCSB and 
superconductivity might appear clearly at /1 ~ 0(A). For the purpose, we intend to make our method parallel with 
the ordinary BCS theory, and supersymmetrize it by using a generalized version of the SUSY Nambu— Jona-Lasinio 
( SNJL ) model [12,13,14,15,16]. We choose the simplest case of the gauge symmetry, namely f7(l), to avoid possible 
complicated situations of breaking patterns of gauge and flavor symmetries. The procedure of this paper will make 
a first step toward the investigation of SUSY superconductivity with more large gauge symmetries. We derive and 
solve a BCS-type self-consistent gap equation of the SUSY superconductivity, and examine the thermodynamics. 
We discuss characteristic aspects of our SUSY BCS theory by comparing with the results of the non-relativistic and 
relativistic BCS theories [1,6,17,18]. To obtain the effective potential at finite temperature, we use the imaginary 
time Matsubara formalism. The Af = 1 SUSY is explicitly broken at finite temperature. The SNJL model was first 
introduced to investigate DCSB in a SUSY field theory [12,13], and it was also applied to the theory of the top 
condensation of the minimal SUSY standard model [8] . The examination on the phenomenon of the chiral symmetry 
breaking of SNJL at finite-temperature and zero-density was given in literature [15]. 
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II. FORMALISM 



The starting point for our investigation is the following Lagrangian of an J\f = 1 SUSY model in four-dimensional 
spacetime: 

C = y'd 2 ^ 2 0(^e y $ + + $ t _e- y $_ + ^$^ t _$ + $_ + ^$^^ + $ + + ^$ t _$ t _$_$_). (1) 

Here, $ + and $_ are chiral matter supcrfields, TV is the number of flavor. G\, G 2 and G3 are coupling constants, 
and they have mass dimension [Mass] -2 . V denotes a real vector multiplet. We consider the gauge degree of freedom 
as U(l)v for the sake of simplicity, and enough for the purpose of this paper. A chemical potential /i is introduced 
by V = 2fid<j°9, as a zeroth-component of vector [9]. Here, the gauge dynamics of U(l) is not considered. In this 
paper, we follow the conventions for metric, gamma matrices, and spinor algebra given in the textbook of Wess and 
Bagger [19]. The third term might give a dynamically generated Dirac mass, while a left-handed and a right-handed 
Majorana masses might be generated by the fourth and fifth terms, respectively. We set aside the question on the 
origin of these nonlinear interaction terms, and introduce them with an assumption of the existence of some effective 
attractive interactions in the system. The special case G2 = G3 ( the left-right symmetric ) will be taken to keep the 
parity symmetry in our model. In fact, the spin-singlet Lorentz-scalar BCS pairing gap is given by a parity-invariant 
combination of a left-handed Majorana mass and a right-handed Majorana mass [17,18]. The original version of the 
SNJL model will be obtained by G2 = G3 = 0. Thus, we call (1) as a generalized-SNJL model. By using (1), we 
examine the dynamical breakings of global U(l)v ( broken by superconductivity ) and U(1)a ( chiral ) symmetries 
at finite temperature and density. 

Our model Lagrangian (1) will be rewritten in the following form through the method of SUSY auxiliary 
fields [13,14,15,16]: 

C = f d 2 6d 2 9( (1 - AW)($t e 2 ^°°* + + Sle- 2 "*' **-) + %-H{H x + £-H%H 2 + ^H^Hs 
J \ G\ G 2 G 3 

+ < y(0)S 1 (£lfi - + 6(0)St(£-Hl - & + &_) 

+S(0)S 2 (^-H 2 - $+$+) +8{6)S\{^Hl - 

+S(0)S 3 (^-H 3 $_d>_) +S (6)Sl(^-Hl d>L4>t_)). (2) 

In order to realize the DCSB in this model, we have introduced a soft SUSY breaking mass A. It was shown 
in the ordinary SNJL model: Because of the nonrenormalization theorem, the chiral symmetry cannot be broken 
dynamically if the theory maintains the Af = 1 SUSY exactly [13]. We speculate it is also the case in (1). Expanding 
the Lagrangian (2) in terms of the component fields, eliminating the auxiliary fields of the chiral multiplets through 
their Euler equations, and keeping only the relevant terms in the leading order in the sense of the large- N expansion, 
we get 



N , . , N , . , N 

' »s 3 



C = -ttM-ttM- Trite 12 



1 <-J2 
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-{d„ - i^)4> ] +{d v + iii5 v0 )te - \ter \ 2 4>lte - Mte 2 \ 2 Ate - A : 
-{d v + iii5 v0 )cj>l{d v - v<W- - \te, \ 2 4>U- - Mte 3 \ 2 <f>U- - A 2 0^- 

-^s^sA-te - l^teA^+te ~ l^te^+te _ Zfiss^sA-te 
-i^ + a v {d v - ifiS I/ o)i/} + - i^)-a v {d v + ifiS u0 )ip- 

+<j> Sl ip+ip- + tfg^+i)- + <j>s 2 ip + ip + + 4>l 2 ^+^+ + tes^P-^P- + tfsji-i)— (3) 

Here, (j} Sj (.7 = 1—3) denote the scalar components of the chiral multiplets Sj ( j = 1 — 3 ). We also assumed that, 
all of cpSj are constant while all of the spinor components of Sj are zero. This assumption is justified in the leading 
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order of the 1/N expansion [13,14,15,16]. Therefore, the partition function is obtained in the following form: 

x (-^l^ 1 | 2 -^l^ 2 | 2 -^|05 3 | 2 ) + 2 i lnDctf! B -zlnDctf> F }, (4) 
where, we have used the following definitions of the fields: 

n-ov —a*), s-r^y s.<»,.*). <« 



Here, \I/ is the four-component Dirac bispinor, S and 5 are the eight-component Nambu notations [17,18,20], and T 
means the transposition. The definitions of the matrices £Ib and f2 F are given as follows: 



/(^+^o)(^+vM-|^ 1 | 2 -4|<^s 2 | 2 -A 2 -2(4,^1+4x^3) 

I -2(4^+4,0*) (d v - itfrtW - itfrt) - \<t> Sl \ 2 - *\<f>s,\ 2 - 



( ^ + 7 V-^S 1 (l + H5)-|4 1 ( 1 -n5) 0s 3 C(l + H5)-4 2 C( 1 -H5) \ ^ 

F V -^ 2 c(i + i 75) + 4 3 c ( 1 - i 75) »p T -7 0T m + ^ 5i (i + h 5 ) + 14, (i - h 5 ) ; ' u 

where, 75 = 7° 7 1 7 2 7 3 , and C = i7 2 7 ° is the charge conjugation matrix. It is a well-known fact that, in the relativistic 
theory of superconductivity, we should take the parity-invariant Lorentz-scalar symmetry for the spin singlet BCS 
pairing state [6,17,18]. In order to take the symmetries of both the Dirac mass and the BCS pairing gap in Oi? 

as Lorentz scalar, we choose 4 X = 0Si an d 4 2 = ^ s ^- ^ 4 2 ^ ^ s 3' ^ e S a P function corresponds to a linear 
combination of the scalar and pseudoscalar pairings. The pseudoscalar pairing cannot realize in the relativistic BCS 
theory [17], and thus we may introduce the constraint 4 2 = &S3 also m the present case. The boson and fermion 
determinants in Eq. (4) will be evaluated in the following forms: 

dem B = (po-Ef(p))(po + Ef(p))(p -E B (p))(p +E B (p)), (8) 
detft F - (p -E F (p)) 2 (p + E F (p)) 2 (p -E F (p)) 2 (p + E F (p)) 2 . (9) 

All of the eigenvalues of S1 F doubly degenerate. These degeneracies relate to the time- reversal invariance of the BCS 
gap function [17]. The bosonic and fcrmionic quasiparticle energy spectra are evaluated as follows: 



E*{p) = \jp 2 + ^ + \4> Sl | 2 + 4\4> S2 1 2 + A 2 - 2vV(p 2 + |0 Sl I 2 + 4|0 S2 | 2 + A 2 ) + 4^ | 2 |</>s 2 | 2 , (10) 
E B {p) = ^p 2 + /i 2 + ^ | 2 + A\4> S2 1 2 + A 2 + 2 v /^ 2 (p 2 + |0 Sl | 2 + 4|^ 2 | 2 + A 2 ) + 4|fe | 2 |^ 2 "F, (11) 
Blip) = ^JWp 2 + \<t>sA 2 -v) 2 + ^s 2 \ 2 , (12) 



E F .{p) = V(vV + S7F + M) 2 + 4|0s 2 | 2 . (13) 

We confirm \<fisi\ is the dynamically generated Dirac mass, while 2|0s 2 | corresponds to the BCS gap function. 

As we mentioned above, we consider the case G2 = G3. The effective action of the leading order in the large- N 
expansion is found to be 

T eff = [ d 4 x(-^^- - 2%^) + 2ilnDctft B -ihiDetQ F . (14) 

Hereafter, we introduce the finite-temperature Matsubara formalism [21]. The Matsubara formalism is obtained by 
the following substitutions in our theory: 

f dpo 1 ■ B ■ F Mr\ 

] M^l^F p Q ^iu n ,iL; n , (15) 
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where (3 = l/fc^T ( ks] the Boltzmann constant, T; temperature ). ks — 1 is taken throughout this paper, ui B and w F 
are the boson and fermion discrete frequencies, respectively. Their definitions are u) B = 2mr/{3 and oj f = (2n + l)7r//3 
( n = 0, ±1, ±2, • • • ). The effective potential ( thermodynamic potential ) becomes 

= %T + + f Wf { E + ip) + EB{P) - E + ip) - 

-| ln(l + e-P E ZW)(l + e-? E -W) + | ln(l - e^^Xl - e"^-^ . (16) 

To obtain the final expression in Eq. (16), the frequency summations were performed. A three-dimensional momentum 
cutoff A was introduced to regularize the integral. The gap equations are derived in the following forms: 







dV, 



eff 



^ - I** I r f (1 - , = ) \ tanh ^ + (1 + ; = ) \ tanh ^) 

G 1 (2tt)3U VFTf^F^f 2 + ^ y/W+W^?' EZ 2 -/ 

+1**1 r^la- ; " 2 + 41052 ' 2 )4coth^ 



fl | M 2 + 4|fe 2 | 2 v 1 coth ^ B ] fiyl 

V vV(P 2 + \^ Sl I 2 + 4|^ 2 1 2 + A 2 ) + 4|0 Sl I 2 I0S 2 | 2 ^ B 2 -/' 



= 



8V. 



eff 



d\4s 2 \ 
4|<fe 2 | 



- 4 '^' / (W{^ tanh f^ + ^ tanh ^} 



G 2 

_gg_/n M 2 + lfeJ 2 )J_ coth ^B 

( 2 +> 3 ^ vV(p 2 + 1^ I 2 + 4|^ 2 1 2 + A 2 ) + 4|^ 5l | 2 |0s 2 F Ef 2 + 
+ (1 + j+^s£ )1 epth^Bl (18) 

At first glance, both of these gap equations seem to include quadratic divergences, and they might be influenced by 
A. The effect of SUSY is parametrized by A. These gap equations correctly give their limiting cases at A — > oo. For 
example, Eq. (18) gives the gap equation of the non-SUSY relativistic superconductivity [17] at A — > oo. Equations 
(17) and (18) can have nontrivial solutions at least at G\,Gi > ( attractive interactions ). The charge density g, 
the conjugate of fi, is found to be 

= 9V ef f 



f d 3 p U- v /p 2 + |0 Sl | 2 , , M+yV + lfeJ 2 , ,P„ F 

I { M 2 E + + w tanh 2 E - 

- ( 1 - P 2 + l^| 2 + 4|^| 2 + A 2 x _^ coth 

V vV(p 2 + 1^ I 2 + 4|0 S2 1 2 + A 2 ) + 4|0 Sl | 2 |0 S2 1 2 J Ef 2 + 

- ( 1 + P 2 + l^| 2 + 4|0s 2 | 2 + A 2 x ^ coth \ 

V vV(P 2 + Ifa I 2 + 4|0 S2 1 2 + A 2 ) + 4|0 Sl | 2 |0 S2 1 2 J E B 2 - 



In the zero-temperature case with \4>s 2 \ = 0> one has 



3 = 2 J (^^-V^ + I^l 2 ) = (20) 
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where, 8(x) is the step function defined as follows: 9(x) = 1 for x > and 6(x) = for x < 0. p F is the Fermi 
momentum. At T = 0, [i coincides with the Fermi energy \/p F + \4>Si I 2 : and it is determined by the charge density 
of fermion. Later, Eq. (18) is solved numerically. We completely neglect the temperature dependence of fi, and use 
it as an external parameter for our numerical calculations. 



III. QUASIPARTICLE EXCITATION ENERGY SPECTRA 

In this section, we examine the quasiparticle excitation energy spectra (10)-(13). Because we consider several 
situations of fi ^ 0, the Bose-Einstein condensation ( BEC ) can take place in our model [9,21]. Depending on the 
model parameters ( A, fi, etc. ) and |0sj and \(f>s 2 \i E+{p) can have a zero-point. In such a case, the logarithmic 

function ln(l — e~ /3E + ) in V e f f will diverge and BEC takes place. A similar discussion on BEC was given in Ref. [9]. 
In this papar, we will not study the physical property of the BEC phase in our model, and only discuss the phase 
boundary in the model-parameter space. This can be done by the examination of the quasiparticle energy spectra. 
We have to find the condition of the realization of BEC before solving the gap equations (17) and (18) to study the 
phases of the DCSB and the superconductivity. 

First, we examine Eqs. (10)-(13) in several limiting cases. At \x = 0, the zero-density case, the energy spectra 
becomes 

E%(p) = VP 2 + I) 2 + A 2 , (21) 

E£(p) = vV + I^P+^l 2 - (22) 
On the other hand, at \4>s 2 \ = one finds the spectra of the DCSB at finite density: 

E B (p) = vV + l0sJ 2 + A2 T/ i, (23) 
E^{p) = vV + l'fePTM, (24) 
while, at |0sj = where only a superconducting gap is generated, the spectra take the following forms: 



E B {p) = v/p 2 +4|^ 2 P + A 2 TAl7 (25) 

E F {p) - v/(IpItm) 2 +4|0 S2 | 2 . (26) 

In Eqs. (23) and (25), E B (p) can become negative with [i > 0. Because of the positiveness of the Bose distribution 
function l/(e^ E + — 1), and /i corresponds to the Fermi energy of the system at zero temperature, 

|05j 2 <^< l0 Sl | 2 + A 2 (27) 

has to be satisfied in Eqs. (23) and (24). If we use [i = \J p F + \^s 1 1 2 , (27) will be rewritten as the condition of the 
external parameters p F and A: 

0<p F <&- (28) 

Thus, A is the upperbound for p F in the SUSY theory. For Eqs. (25) and (26), 

0<A 1 2 =p|<4|^ 2 | 2 + A 2 (29) 

has to be satisfied. From (23) one finds that, if dynamically generated |0Si I obtained as a solution of the gap equation 
(17) satisfies /i 2 = |05- x | 2 + A 2 , the BEC takes place at the mode of p = 0. 

Figure 1 shows a typical case of the excitation energy spectra of boson and fermion qu asiparticlcs under the 
superconducting state with a non- vanishing chiral mass. E F (p) has a minimum at p F = -\//i 2 — \<t>Si I 2 where the 
excitation energy gap 4|</>,s 2 | locates. The appearance of the branches E B {p) and E B (p) is the new phenomenon of 
our SUSY theory compared with the no-SUSY relativistic BCS superconductivity [6,17,18]. The Bose branch E B has 
an energy gap at \p\ = 0. All of the spectra become parallel with the light cone at |p| — > oo. Let us examine the 
situation E B (p) = in detail. From Eq. (10), one obtains the solutions of the equation E B {p) = as 



\p\ - vV - I^Sil 2 - 4|<As 2 | 2 - A 2 ± 4|0 Sl pal . (30) 

From the positiveness of the Bose distribution function, we have to avoid the situation where Ef(p) = has two 
solutions, because any |p| inside these two solutions makes E B (p) a complex number. Therefore, we consider the 
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following three cases in Eq. (30): (I) |0sj = and \<f>s 2 \ ( superconducting state without DCSB ), (II) \(f>s 2 \ = 
and \(psi\ ( DCSB without superconductivity ), (III) |0 Sl | = \<fis 2 \ = 0- I n the case (I), Eq. (30) becomes 
\p\ = \/ — 4|(/>s 2 | 2 — A 2 . From this solution, one finds the condition fi 2 > 4|0s 2 | 2 + A 2 nas to be satisfied. Thus, 
taking into account (29), we get the critical SUSY soft mass of the superconductivity: 

A£< J > = vV-4|fe 2 | 2 ~PF. (31) 

Here, we used fj, = pf- When A = Acr^ in the case (I), Ef(p) has one zero-point at \p\ = 0, and it behaves 
at \p\ — + as Ef(p) w p 2 /2\J (Ac£^) 2 + 4|0s 2 | 2 - Thus we conclude that, in the case (I), the superconductivity 

can coexist with the BEC at A = A^ J \ while only the superconductivity can take place at A > A s C r^ ■ In other 
words, the phase boundary given by (31) depends on a self-consistently determined \<fis 2 \- I n the case (II), we find 
\p\ = \/ n 2 — \4>S! I 2 — A 2 as one solution of Ef(p) = 0. From this solution, one has [i 2 > |0sj 2 + A 2 , and with the 
condition (27), we obtain the critical SUSY soft mass of the superconductivity: 

A^ (77) - V^-\^ Sl \ 2 =PF. (32) 



Here, we used fj, = \fp 2 F + \^Si | 2 . When we choose the model parameters suitably, especially to satisfy A > Acr 
with \4>s 2 \ 7^ 0; Ef(p) takes a positive value for any |p| as shown in Fig. 1. On the other hand, the BEC in the boson 

sector takes place at A = A s c f H \ and E+(p) behaves as Ef(p) rts p 2 /2\J (A s c C r Il) ) 2 + |0sj 2 at \p\ — » ( there is no 
energy gap of the branch Ef ). In this situation, \(f>s 2 \ = has to be satisfied, thus the superconductivity cannot 
coexist with the BEC. We find that the chiral symmetry ( chiral mass {(fisA ) affects the physical situation of the 
phase boundary between the superconductivity and the BEC. These results are summarized in table I. 



IV. DYNAMICAL CHIRAL SYMMETRY BREAKING 



Now, we arrive at the stage to study the DCSB without superconductivity. At T — and \(j>s 2 \ = 0, the gap 
equation (17) becomes 

Gi [ A 2 , ( 0(p-Pf) 1 



1 _ fca / 2 d ( V{P-Pf) i \ 

2 W P p Vp 2 + I^I 2 Vp 2 + ttsA 2 + A 2 ) 

\f\4>Si\ + A 2 

PF + VPf + IfaA 2 



-Pf^JpI + \4> Sl I 2 + \4> Sl 1 In l r ^-j j • (33) 

With taking the limit A — > oo, Eq. (33) gives the gap equation of the DCSB of the ordinary ( no-SUSY ) Nambu— Jona- 
Lasinio model at finite density [5] , while at pp — it will give the gap equation of the zero-density case of the SN JL 
model [13,15]. This equation determines the phase diagram in the space of parameters GiA 2 , A/A and pp/A. The 
examination of the critical coupling might be the easiest method to see the relation of several possible phases. The 
determination equation for the critical coupling (Gi) cr is obtained as follows: 



4vr 2 

(Gi) cr A 2 = ■ . (34) 

A 2 A 2 i+yi+A'/A' _p% 

F t F m A/A A 7 







Figure 2 shows (Gi) cr A 2 as a function of A/A. Equation (34) gives (Gi) cr A 2 of the no-SUSY case as 
lim A / A ^ 00 (Gi) cr A 2 = 47r 2 /(l — pp/A 2 ). The denominator of Eq. (34) includes the finite-density effect on (Gi) cr . 
(Gi) cr of pf is larger than that of pp = 0. The divergence of (Gi) cr A 2 at a non-zero value of A/A ( depends on 
a numerical value of pf ) indicates the existence of the critical soft mass A^ sb for the DCSB of this model [13,15]. 
From (28), we know pf < A has to be satisfied in Fig. 2: When pp — A, the BEC takes place at the mode of p = 0. 
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V. SUPERCONDUCTIVITY 



Now, we consider the following situation: After the Dirac mass \<j>s 1 1 dynamically generated and the chiral symmetry 
was broken, we assume the superconducting instability occurs in the system. We solve the gap equation (18) under 
this situation, and restrict ourselves to examining the superconductivity. Thus, |05' 1 1 is treated as a model parameter. 
Usually in various superconductors we know, the BCS gap function 2|0s 2 | is much smaller than fj, — l^gj [18]. We 
assume it is also the case in our numerical calculation for solving Eq. (18), and choose the model parameters G2, A, fi, 
A and l^sA to satisfy 2|0s 2 |/( / u — l<Asil) ^ 1- Because 2|0s 2 | is quite sensitive to G2 and A even in our SUSY theory, 
we carefully choose numerical values of them ( a kind of fine-tuning ) to get a physically reasonable solution. In this 
paper, we set aside the question on the origin of the interactions given from the nonlinear terms of <I>±. We treat 
the cutoff in a general way. As discussed in Sec. Ill, the superconductivity never occurs at pf — A with |(/>sj 7^ 0. 
We also have to take into account the condition (28) to choose a numerical value of A; when A < Acr'^ = Pf is 
satisfied, the superconductivity cannot realize. 

Figure 3 shows the gap function 2|</>s 2 | at T — as a function of the coupling constant Gi. Equation (18) was 
solved under the three examples: A = 2, A = 10 and A = 100 with the energy unit |</>sj = 1- The result of 
the no-SUSY case ( obtained by dropping the contribution coming from the bosonic part in Eq.(18) ), namely the 
ordinary relativistic superconductivity, is also given in this figure. All of the examples show qualitatively similar 
dependences on G2, and they reflect the nonperturbative effect. When the soft mass becomes large, the dependence 
on G2 becomes strong, however, the dependence is qualitatively unchanged whether a numerical value of A is larger 
or smaller than A. The result in Fig. 3 also indicates the absence of the critical coupling of the superconductivity 
similar to the ordinary BCS theory [1,6,17,18]: Any attractive interaction gives a superconducting instability also in 
the SUSY case. In Eq.(18), the contribution coming from the vicinity of the Fermi energy dominantly determines 
the solution |^>s 2 |- This contribution can be made arbitrarily large by changing \4>s 2 \- The integrals given from the 
momentum- integration of the functions of the branches E F , E? and E® are almost constants under the variation of 
|(/f>s 2 |. Thus, Eq. (18) at T = always has one nontrivial solution with any attractive interaction. We also find that, 
Eq. (18) never has two solutions in our numerical calculations. In other words, the effective potential (16) only has 
one minimum with respect to the variation of \</>s 2 \ under the superconducting phase. We confirmed this fact by our 
numerical calculation of the effective potential. 

In Fig. 4, we show the gap function 2|0g 2 | at T = as a function of the three-dimensional momentum cutoff A. In 
the case A = 2, the divergence of the gap function becomes slower than other examples because of the effect of SUSY. 
In this figure, the gap function at A = 2 diverges almost linearly. In the cases of A = 10, A = 100 and no-SUSY, the 
gap function diverges almost quadratically at A — > 00. 

Figure 5 gives 2|0s 2 | as a function of temperature T. 2\<j)s 2 \ continuously vanishes at T — > T c ( T c ; the critical 
temperature ), clearly shows the character of second-order phase transition. In all of the examples shown in this figure, 
the BCS universal constant 2\<j) S2 (T = 0)|/T C = 1.76 is satisfied [1,18]. This fact means that, both 2\<ps 2 { T = 0)| and 
T c depend on G2, A, A and /1, while their ratio is independent on them. In the non-relativistic and relativistic BCS 
theories, any spin-singlet ( or, Lorentz-scalar ) BCS gap function satisfies the universal constant [1,18]. We find it is 
also satisfied in our SUSY BCS theory. 

We discuss the thermodynamic property of the SUSY superconductivity. With taking into account dV e f / / 'd\(j>s 2 \ = 
0, the entropy S is obtained as follows: 



q dV eff 

• f A d3p ( An ( l + e ~ 0E +)( l + e ~ l3E -) P E + , & E - , $ E + , p eB 



(2tt) 3 I (l-e-^Xl-e-^-) e^ + l eP E - + 1 e^-l e? E - 



-}■ 05) 



This is the entropy of the ideal gas of quasiparticles of the branches E F , E F , Ef and E 1 ^. At \4>s 2 \ — * 0, (35) gives 
the entropy of the normal state. By the result shown in Fig. 5, we find S is continuous at the phase transition: It is 
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a second-order phase transition. The heat capacity C becomes 




Let us recall the energy spectra shown in Fig. 1. Because the density of states diverges around the energy gap of 
[18], the contribution of the excitation of E+ dominates the integrals of S and C, and other branches give almost 
no contribution. Both the entropy S and the heat capacity C of the system are determined by the thermal excitations 
of the quasiparticles of the branch E^. Because the temperature dependence of 2\(j>s 2 | is the same with the well-known 
BCS result, we conclude that the temperature dependences of S and C are qualitatively the same with the well-known 
results of the ordinary BCS theory [1,18]: The heat capacity C of Eq. (36) becomes exponentially small in the limit 
T — > 0, while it behaves as C oc T ( the Fermi liquid behavior ) at T > T c . 



VI. SUMMARY AND DISCUSSION 



In summary, by examining the quasiparticle excitation energy spectra and the gap equations, we have discussed 
the DCSB and the superconductivity in the generalized SNJL model. We have found the finite-density effect in the 
critical coupling of the DCSB. The effects of the bosonic part in the gap equation (18) in the SUSY BCS theory 
have been reviewed in detail, while we have not found any SUSY effect in the thermodynamic character of the SUSY 
superconductivity. These results have been understood by the energy spectra (10)-(13). We have revealed that, the 
superconductivity shows the BCS character even if \x is close to A. There is no critical coupling also in the SUSY BCS 
case, while we have found the critical soft mass Ac£^: SUSY protects not only the chiral invariance from DCSB but 
also the dynamical breaking of gauge symmetry at finite density. If the soft mass A is slightly larger than A^ 11 ^, 
the SUSY effect becomes significant in the magnitude of the gap function 2|(/>g 2 |. 

Finally, we would like to make some comments on several issues and possible extensions of this work. It is interesting 
for us to examine the collective modes and its excitations in both the DCSB and the superconductivity. An interaction 
between bosons might alter the excitation energy spectra of the boson sector, as discussed in the Bogoliubov theory 
of superfluidity [22]. 

The existence of the upperbound of fj, for the DCSB and superconductivity is a remarkable fact, when we consider 
phenomenological aspects or cosmological problems. Today, the superpartners of the known elementary particles are 
supposed to exist in the TeV energy scale. From this point of view, any matters of massive Dirac particles cannot take 
the chemical potential over the TeV region without destroying the superconductivity. In this paper, we choosed the 
gauge freedom broken by the superconductivity as Z7(l) for the sake of simplicity. We regard our result is the starting 
point to extend the theory to the more general SU(N C ) case. From the phenomenologically possible energy scale 
of A ( >TeV ), our theory is relatively closer to the electroweak theory, especially the top condensation model [8], 
or the technicolor theory [7] than to SU(3 C ) QCD. For example, the SU(3 C ) gauge interaction itself is the origin of 
the attractive interaction for CSC. Therefore, if one uses an NJL-type model to describe CSC, the model parameters 
should be chosen from the consideration of the QCD gauge interaction and hadron phenomenology. A coupling 
constant for an NJL-type model should take a numerical value of C(GeV" 2 ), while a cutoff will become O(GeV) [5,6]. 
At the order of the energy scale of the cutoff, there is no SUSY effect in CSC in the sense of the context of this paper. 

The applications of our method to an investigation of the SUSY superconductivity in (2+l)-dimensional case [23] can 
easily be done. The gauge dynamics was completely neglected in this paper. The gauged-NJL model was extensively 
studied in the context of the strong-coupling QED [4]. Effects of gauge fields on the DCSB and superconductivity 
can be examined by our model at finite temperature and density. The extensions to several gauge groups of SUSY 
grand unified theories ( SGUTs ) are possible for our theory [24] . Recently, the phase structures of the ordinary N JL 
model with an external electromagnetic field, both the cases in curved spacetime [25] or at finite temperature and 
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density [26], are studied. The DCSB of the SU(N C ) SNJL model at zero density was examined in several curved 
spacetimc to describe the physics of the early universe. The extensions of our theory of this paper to these several 
external conditions are interesting [27]. 

In this paper, we have given the possibility of the calculations of several thermodynamic quantities ( thermodynamic 
potential, entropy, heat capacity ) in a SUSY condensed matter. It is one of the advantages of our method. From our 
method, the Ginzburg-Landau theory for the DCSB and superconductivity can be constructed and applied to study 
vortex in the SUSY case. Our theory provides a way to study the SUSY condensed matter physics. The extension of 
our theory to the case of a vector-like SU(N c )-gauge model with A^-flavor will be published elsewhere [28]. 
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FIG. 1: The excitation energy spectra of quasiparticles under the superconducting state with a finite chiral mass. We choose 
[i = 2, \<t>s 2 \ =0.02, A = 2 and \<j> Sl \ = 1. 

FIG. 2: The critical coupling (Gi) C rA 2 of the DCSB, shown as a function of A/A. 

FIG. 3: The gap function 2\(j>s 2 1 at T = 0, shown as a function of G 2 - We set n = 2, A = 5 and \(f>s 1 \ — 1. 

FIG. 4: The gap function 2|0g 2 1 at T = 0, given as a function of A. We set /u = 2, G2 = 0.5 and |<^>s 1 = 1. 

FIG. 5: The gap function 2\<f>s 2 \ given as a function of T. The model parameters are set as fi — 2, A = 5, G2 = 0.5 and 
\4>s x \ = 1- The curves, taken from the uppermost one, correspond to the values of A = 100, A = 10, and A = 2 in the unit 
1**1 = 1- 



TABLE I: Relations of several orders. 
SUSY breaking mass BEC DCSB CSC 



A = \J p? — |0i| 2 takes place possible no solution 
A = \//i 2 — 4|02p takes place no solution possible 



Erratum: Dynamical Chiral Symmetry Breaking and Superconductivity in the 
Supersymmetric Nambu— Jona-Lasinio Model at finite Temperature and Density 



Tadafumi Ohsaku 
Institute for Theoretical Physics, University of Cologne, Germany 
(Dated: April 20, 2008) 

Some equations of the paper, T. Ohsaku, Dynamical Chiral Symmetry Breaking and Superconduc- 
tivity in the Supersymmetric Nambu— Jona-Lasinio Model at finite Temperature and Density, Phys. 
Lett. B634 (2006) 285-294, must be corrected. Fortunately, the numerical results of the original 
version of this paper are still valid for our consideration on the supersymmetric BCS superconduc- 
tivity. 



PACS numbers: 11.30.Pb, 11.30.Rd, 74.20.Fg 

bo 

| First, several equations we must correct will be given bellow. The matrix Qp of Eq. (7) in the original paper must 
£vq . be replaced by 

F V -0&Cr(l + iT6)-4 8 C7(l-i<ys) if - 7 0T M + \<t> Sl (1 + Hz) + |4 X C 1 " *7 5 ) J ' 

Here, the sign of the BCS-gap functions <frs 2 an d 4>s 3 of the block-off-diagonal part of is different from that of the 
original version. We derived the fermion quasiparticle excitation spectra and the fermionic contribution inside the 
self-consistent gap equations correctly in the original version, while we must correct several formulae of the bosonic 
sector. All of the places we mush correct are coming from this modification of signs of the gap function inside ftp. 



> 



C3 



In order to take the symmetry of the BCS pairing gap in £lp as a Lorentz scalar, we should choose now, 



ON . 

( We have chosen $ s — 4>s 3 in the original version. ) Under this choice, the bosonic quasiparticle energy spectra of 
Eqs. (10) and (11) of the original version will be replaced by 



E*(p) = VV + l0s 1 | 2 +4| ( fe 2 | 2 + A 2 -/i, 
^ EB (P) = V / P 2 + l0s 1 | 2 +4|fe 2 | 2 + A2 + / i. 

^ '" Note that E± depend on p 2 as monotonically incleasing functions. Therefore, the gap equations Eqs. (17) and (18) 
C"| , of the original version are replaced as 



= 



dV f 



X 

2J ^ " \fa I / A { (1 - = ) -4 tanh £f£ + (1 + = ) X tanh V) 

t £l 1 [ coth p_ e b coth p_ e b\ 



o = 



+1**1 

BVeff 

d\<Ps 2 \ 
G 



2 



4| ^'/ (2^{^ tanh ^ + ^ tanh ^l 
+4 U S J f f* - 1 (coth P-E*+ coth ?-E B ), 

respectively. Now, the charge density g given in Eq. (19) in the original paper becomes 

= dV e ff 

jPp ( tanh P_ EF + m + v^TTCT tanh V + coth V - coth V 

(2tt) 3 I £f 2+ £ F 2 " 2 + 2 " 



2 



The energy spectra E± at \i = given in Eq. (21) of the original paper must be corrected into 

Ei( P ) = v / P 2 + l^5 1 | 2 + 4|0s 2 | + A 2 . 
The solution of the equation Ef (p) = given in Eq. (30) of the original paper will be replaced by 

\P\ = vV-lfeP-^P-A 2 - 

We give the condition on the Bosc-Einstcin condensation ( BEC ) of the scalar sector by the expression of Ef given 
above. For the positiveness of the Bose distribution function l/(e f3E + — 1), 

Vl0 Sl | 2 +4|</> S2 | 2 + A 2 >M>0 

must be hold. The BEC takes place when /j, = y/\<j)Si | 2 + 4|(/>s 2 1 2 + A 2 is satisfied. This relation includes (27) and 
(29) of the original version. Therefore, the existence of the upper bound of fj, ( ~ 0(A) ) is true. While, [i is the 
Fermi energy at the zero-temperature case, which we defined in the original version of this paper as follows: 

0< \<p Sl \ <M- 

Hence both of these inequality relations must be satisfied, as we revealed in the original version. 

Finally, we wish to mention the numerical results given in the Sec. V of the original version. ( We do not have 
to correct the numerical result of Sec. IV in the original paper. ) The results are still valid. In Ref. [1], the author 
discussed a similar system, a generalized Supersymmetric Nambu— Jona-Lasinio model with non-Abelian gauge/flavor 
degrees of freedom. The numerical results of Ref. [1] ( solving the gap equations ) were performed under the following 
two cases: (i) dynamical chiral symmetry breaking without color superconductivity, (ii) color superconductivity with 
a vanishing chiral mass. Hence, the numerical results of Ref. [1] did not consider the situation where two orderings 
coexists ( i.e. the case |0sj ^ 0, \<fis 2 \ )• I n f act 7 errors inside the gap equations of the original version of 
this paper arises when we consider the case of coexistence of two orderings. Furthermore, the numerical results of 
the original version of this paper are qualitatively the same with that of Ref. [1] . The reason of the validity of the 
numerical results of the original version is that usually a contribution of the boson sector in a gap function is small 
compared with that of the fermion sector. ( We confirm that numerical errors in Fig. 1 are 0(1O~ 3 )-0(1O -4 ), due to 
our choice 1 0s 2 1 / 1 4>s 1 \ — 0.02. In fact, all of branches E±, E± have gaps, and several thermodynamic quantities are 
dominated by the contribution coming from the branch E+ because the density of states of E+ diverges at |fc| = k F . 
) Thus, numerical errors of the gap functions in the original paper are quite ( invisibly ) small. Hence, we conclude 
that the BCS character of the supersymmetric superconductivity we obtained in the original version is also completely 
valid. 



[1] T. Ohsaku, Dynamical Chiral Symmetry Breaking, Color Superconductivity, and Bose-Einstein Condensation in an 
SU(N C ) x U(Nf)L x U '(N f) R-invariant Supersymmetric Nambu— Jona-Lasinio Model at finite Temperature and Density, 
hep-ph/0601009, submitted for publication. 



